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We present a scheme to conditionally generate an arbitrary superposition of a pair of coherent 
states from a squeezed vacuum by means of the modified photon subtraction where a coherent state 
ancilla and two on/off type detectors are used. We show that, even including realistic imperfections 
of the detectors, our scheme can generate a target state with a high fidelity. The amplitude of the 
generated states can be amplified by conditional homodyne detections. 
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I. INTRODUCTION 

Conditional quantum operation based on photon de- 
tection plays an important role in recent optical quan- 
tum information processing. Particularly, in the contin- 
uous variable regime, it is the only available tool with 
current technology to generate non- Gaussian states from 
Gaussian ones. A typical example is the 'photon sub- 
traction' operation, in which a nonclassical input state 
is split by a highly transmissive beamsplitter (BS) and 
the reflected state is measured by a photon number re- 
solving detector (PNRD). Selecting the event that the 
detector observes photons, one obtains a non-Gaussian 
transformation from the input to output quantum state. 

This type of conditional operation was formulated in 
Dakna et al. [2] then showed that applying it to a 
squeezed vacuum, one can generate a non-Gaussian state 
which is close to the superposition of coherent states with 
plus or minus phase 

|C±(a)) = -^(|a)±|-a)), (1) 

with very high fidelity, where \a) is a coherent state with 
the amplitude of a and A/± is the normalization factor. 

Recently, such state has been experimentally gener- 
ated by a single photon subtraction from pulsed [3|, [J] 
and CW [sl, [oj squeezed vacua. In these experiments, 
since a reflected beam includes sufficiently small average 
number of photons (n ^ 1), single photon detection was 
approximately realized by use of an avalanche photodi- 
ode (APD) which is often called as 'on/off' type detector 
since it discriminates only a presence of photons instead 
of resolving photon numbers. 

The progress of these experiments promises the realiza- 
tions of more complicated applications of photon subtrac- 
tion proposed so far, including the improvement of quan- 
tum teleportation [7, 8, 9] and entanglement-assisted 
coding [10], entanglement distillation [II], loophole free 
tests of Bell's inequalities [H, [isl, and 02tical quan- 
tum computations in quadrature basis [lEj, or super- 
posed coherent state basis [17, 18]. In the last applica- 
tion, \C±{a)) with appropriate a is required as an an- 
cillary state. To prepare such ancillae, the method to 



conditionally amplify a with on/off detectors has been 
proposed [T3, [l^. Fiurasek et al. [2l| also recently 
showed that one can arbitrarily synthesize a single-mode 
quantum state up to the A/'-photon eigenstate by concate- 
nating squee zing operations and N times single photon 
subtractions [2l|. 

In this paper, we propose a method to conditionally 
generate the state in which two coherent states are su- 
perposed with arbitrary ratio and phase, c+|a) H-c_ | — a). 
This is accomplished by a simple modification of the 
scheme proposed by Dakna et al. (DAGKW) [2] . We first 
discuss an ideal scheme using two PNRDs and a qubit 
ancilla, which produces a superposition of the one- and 
two-photon subtracted states. We show that such state 
fairly well approximates the target state c+ |a) +c_ | — a). 
We next present a more practical scheme where PNRDs 
and a qubit ancilla are replaced by the on/off detectors 
and a coherent state ancilla. Even including practical 
imperfections of the detectors, it can generate the target 
state with a high fidelity. 

Our scheme should be compared with the one by 
Fiurasek et al. [21], which requires A^ detectors to syn- 
thesize a state consisting of the number states up to lA"). 
Ours, on the other hand, uses only two detectors to syn- 
thesize a fully continuous variable state while, in return, 
the class of states to be generated is restricted. 

Finally, we show that our scheme is useful to simplify 
the setup of the conditional amplification of the superpo- 
sitions of coherent states originally proposed in [1^, [20] • 

The paper is organized as follows. In Sec. [Ill we dis- 
cuss an ideal setup with PNRDs and a qubit ancilla and 
how the scheme's parameters are optimized to generate 
desired superposed states. In Sec. IIIIl a practical scheme 
using on/off detectors and a coherent state ancilla is 
shown and its experimental feasibility is numerically ex- 
amined. An application of our scheme to the conditional 
amplification of superposed coherent states is shown in 
Sec. [IV] and Sec. [Vj concludes the paper. 
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FIG. 1: (a) Generation of a plus- or minus-superposition of 
coherent states via photon subtraction operation with a pho- 
ton number resolving detector (PNRD). (b) Generation of an 
arbitrary superposition of coherent states with PNRDs and a 
qubit ancilla. 
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FIG. 2: Fidelities between the photon subtracted states and 
the ideal superposition of coherent states, (a) (b) 
KC-(a)|^l)|^and (c) |(C+(a)|^2)r- 



II. GENERATION OF AN ARBITRARILY 
SUPERPOSITION OF COHERENT STATES 

Figure [ija) illustrates the DAGKW photon subtrac- 
tion scheme [2]. A squeezed vacuum with the squeezing 
parameter r is mixed with a vacuum by a highly trans- 
missive BS and the reflected part of the state is detected 
by a PNRD. When the reflected part is projected onto 
the photon number eigenstate |m) (m > 0), the state re- 
mained in the transmitted mode is reduced to be the m 
photon subtracted squeezed vacuum state that can be de- 
scribed by a minus- or plus-superposition of two distinct 
states as 



ivi/™) = A(ivi/L+)) + (-iri*L-))), 



(2) 



where A is the normalization factor [22.]. It was shown 
that, with an appropriate input squeezed vacuum, the 
states 1^^^) are very close to the coherent states | ± 
am) and thus the states |^rn) are also very close to a 
superposition of coherent states. 

Let us extend the above scheme as illustrated in 
Fig.fljb). Let the upper BS has the power transmittance 
T ^ 1 and the lower be a balanced BS. The reflected part 
of the state is mixed with the auxiliary state 6o|0) + 
and then each port is incident into a PNRD. After some 
calculations, one finds that if the measurement outcome 
of the two detectors is (2, 0) or (0, 2), the reflected part is 
effectively projected onto ^bl/V2\l) +6o/V2|2) and the 
transmitted state conditioned on either of these outcomes 
has the form 



l^out) = ail^i) +^21^2), 



(3) 



where ai and a2 are the functions of 6o, 6i, and T. Since 
l^rn) can be regarded as a superposition specified in 
Eq. ([1]), the state |^out) is also expected to be a super- 
postion of \ ±a) with the controlled ratio and phase by 
choosing ancilla parameters bo and hi appropriately. 

Now let us see the state in Eq. ([3j) more carefully. 
To show l^out) to be a superposition of two (classical) 



macroscopically distinct states, one has to find the fol- 
lowing decompositions. 



1*2) = + 



(4) 
(5) 



in which \<p±) are close enough to the coherent states 
I ci and C2 are the normalization factors satisfying 

\(j)±) = C2|^2) Note that the decomposition 

described in Eq. ([2]) [2] is not optimal in our purpose 
since and |^2) do not share the common decom- 

posed components. 

The optimal \(j)±) to maximize the fidelity |(^±| 
can be derived from the exact expression of |^rn) [22] and 
Eqs. (HI) and ^ as 
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where A = tanhr is the squeezing parameter and the 
amplitude of the corresponding coherent states is given 

by 
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Then we have quasi-coherent states 

(1 - A2T2)3/4 
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and the fidelity between Eqs. ( 
= Vl- A2T2(1 + AT)(1 



and (|9]) is given by 



2AT) exp 



3AT 



1 + AT 



(10) 

which is plotted in Fig. ([2]) by the black line (line (a)). 
For a < 1, more than 0.99 fidelity is achieved. The va- 
lidity of this optimization is also confirmed by looking 
at the fidelities |(^i|C_(a))p and |(^2|C+(a))p for the 
same a. These are plotted in the same figure by the red 
(line (b)) and green (line (c)) lines, respectively. 
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FIG. 3: A schematic of the generation of an aribitrary super- 
position of coherent states with on/off detectors and a dis- 
placement operation. 



III. PRACTICAL SETUP WITH ON/OFF 
DETECTORS 



Preparing PNRDs and a photon number qubit an- 
cilla is still somehow challenging with current technology. 
In this section, we show a modified and more practical 
scheme in which PNRDs and a qubit ancillary state are 
replaced with on/off detectors and a coherent state, re- 
spectively. 

The modified scheme is depicted in Fig. [3l The re- 
fiected state from the first BS with the transmittance T 
is further split by the second balanced BS. One beam is 
directly measured by an on/off detector (mode B) and 
the other is first shifted by the displacement operator 
D{P) = exp[pd^ — p*a] and then measured by another 
on/off detector (mode C). It is well known that a dis- 
placement operation is realized by interfering the signal 
with an auxiliary coherent state \/3/\^l — Tjj) by a BS 
with the transmittance Tjj. In the limit of Tjj 1, this 
operation is exactly the same as D{P). The output state 
is conditionally selected only when both detectors are si- 
multaneously clicked by photons. The photons detected 
in mode B always come from the squeezed vacuum while, 
in mode C, the photons from the squeezed vacuum are 
interfered by the displacement. This quantum interfer- 
ence and the on/off detection realizes a projection onto 
a superposition of different photon number states. 

The positive operator- valued measure (POVM) for 
on/off detectors is described by {Iloff, Hon} where Iloff = 
|0)(0| and lion = I — Hoff and / is an identity operator. 
Similarly, when a displacement operation D{/3) is placed 
before the detector, as in mode C, the total on/off POVM 
is expressed as 

ilosiP) = DHm){O\D{0) = \-P){-P\, (11) 

no„(/3) = D\p){i-\0){0\)D{p) 

= (12) 

The average photon number refiected to mode B from 
the initial squeezed vacuum is given by (1 — T) sinh^ r. 
For moderate squeezing, this is sufficiently small to as- 
sume that the reffected beam in mode B contains maxi- 
mally one photon and to ignore the more than one pho- 
ton eigenspace at the measurement process (e.g. (1 — 



T)sinhV - 0.005 for r = 0.3 and T = 0.95). When 
in Eqs. (pT]) and ([12]) is also sufficiently small such 
that one can approximate as | — /3) ^ |0) — Iloff (/?) 
in mode C acts as the projection onto |0) — and 
i^on{/3) as the projection onto the orthogonal superposi- 
tion /3*|0) + As a consequence, when both detectors 
are clicked, the reflected part of the state is projected 
onto 



\Bl^,\l)B{f3*\0)c 



c\v\n,,^\^/B {P |u;c + |l)c) 
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/3*\1)b + \2) 



B, 



(13) 



where the normalization factors and global phases are 
omitted and Bt = exp[0{a^b — ab^)] is a BS operator 
with cos^ 

On the other hand, the state after the flrst BS is de- 
scribed as 



BT[s{r)\0)A)\0)B = 



Po|*o)a|0)b + 
^21*2)412)3 



Pi|*i)a|1)b 
(14) 



where S{r) = exp[r/2{a^ — a^^)] is a squeezing operator 
and Vm is the probability to observe m photons in modes 
B andC 0, 
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(15) 



From Eqs. (p!3|) and (p!4|) . we have the conditional output 
for the simultaneous click as 



|*out) OC/3v/:Pi|*i) + \/^2|*2). 



(16) 



Consequently, for the generation of the superposition 
state c+|^+) + c_|^_), the optimal displacement 13 is 
derived from Eqs. (|4l-[7]), (p!5|) . and (p!6|) and given as 



c+ 



3A(i-r) V'^ 
2(1 - A2r2) ) 



(17) 
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which is vahd under the condition of 



< 1, i.e. 



2(1 - A^T^) 
3A(1 -T) 



(18) 



Note that Eq. ([T7j) is almost optimal for arbitrary P al- 
though this condition will be broken when c+ + c_ ~ 
i.e. one wants to generate |^out) ^ l^i)- For large 
Eq. (p!2|) up to one photon state is given by 



non(/3) - (1 - e-l'" )/ + e-l''l (/3*|0) + |1))(/?{0| + {1|). 

(19) 

Although it makes the output as a mixed state of 

Pout = (1 - e-l^l')|M/i)(*i| + e-l'^l'|vl/out)(*out|, (20) 

this is clearly a negligible error since |(^out|^i)P expo- 
nentially approaches to unit. 

In the rest of this section, we numerically examine the 
conditional outputs under realistic conditions. In prac- 
tice, there is always finite probability to detect more than 
one photon at each detector. Moreover, the detectors 
themselves have finite imperfections. The POVM for an 
imperfect on/off detector with the displacement opera- 
tion D{a) is given by 

oo 

flos{a,v,'y) = e-^'Y^il -vrD\a)\m){m\D{a), 



m=0 



non(a,?7, zv) = /-noff(a), 



(21) 
(22) 



where 77 and u are the quantum efficiency and dark count 
of the detector, respectively. 

This kind of detectors makes the output an unwanted 
mixed state. To derive photon subtracted states under 
these conditions, it is useful to use the characteristic func- 
tions to describe the states and POVMs [H, HI, |25|. 
Since the input squeezed vacuum is a Gaussian state, 
its characteristic function can be described as 



Xsv(^) = exp 



1 



(23) 



where co = iu^v)^ is a two dimensional vector and Fsv 
is the covariance matrix for the squeezed vacuum 



■ sv 









,-2r 



(24) 



The mixing of a squeezed vacuum and a vacuum by a BS 
is described by a linear transformation 

ESV © Evac ^ Sls{T)TsY © Evac^B5(T), (25) 

where Fyac = I is the covariance matrix for the vacuum 
state and Sbs{T) is the 4x4 matrix 



Sbs{T) = 



Vti 



Ti Vti 



(26) 
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FIG. 4: The Wigner functions of the superpositions of co- 
herent states (c+|(/)+) + c-\(j)-)) / J\f generated by the on/off 
detector photon subtractions illustrated in Fig. (|3]) for (a), 
(b) {c+,c-} = {l,z}, (c) {c+,c-} = {3,-1}, and (d) 
{c+,c_} = {1,0} where (a) r = 0.3, T = 0.999, a = 0.97 
with the ideal detectors of 77 = 1, and = 0, and (b), (c), (d) 
r = 0.3, T — 0.95, a — 0.95, with the practical imperfection 
of detectors, 77 = 0.1, and v — 10~^. The fidelities between 
the plotted state and the ideal state (c+|a) + c-| — a))/J\f 
are (a) F = 0.993, (b) F = 0.952, (c) F = 0.978, and (d) 
F 0.994. 



Then the covariance matrix after the two BSs in Fig. [3] 
is given by 



{I®Sls{l/2)){Sls{T)®msw 
x{Sbs{T)®I){I®Sbs{1/2)). 



(27) 



Let the characteristic function corresponding to F be 
Xin(^A, ^B, ^c)- The output characteristic function con- 
ditioned on the simultaneous click in both two detectors 
is then given by 

/CO fOO 
/ dc^B dcc;c Xin (^ a , , ) 
-00 J — 00 

XXon(-^B,-Cc;c), (28) 

where 

Xon(^B,^c) = Xon(0,CJB)Xon(-/3*,C^c), (29) 

and Xon(<^,^) corresponds to the characteristic function 
for the POVM defined in Eq. ([22]). Finally, the Wigner 
function for the output is given by the Fourier transform 
of the characteristic function as 

1 

^out(z) = J dcJXout(^)exp [-zcj^z] , (30) 



where z = (x^p)^ . 

Figure m^a) shows Wout(z) corresponding to the state 
with nearly ideal parameters, T = 0.999, unit 
quantum efficiency, and zero dark counts. The squeezing 
parameter r = 0.3 corresponds to the coherent ampli- 
tude of a = 0.97. The fidelity between this output and 
\a) -\- i\ — a) is F = 0.993. Note that this is the same 
as the state generated from a coherent state by strong 
Kerr nonlinear evolution [26]. Figures [Hb)-(d) plot the 
Wigner functions with realistic detectors (77 = 0.1 and 
u = 10~^) for different {c+,c_}. Even with such imper- 
fect detectors, high fidelities (> 0.95) could be achieved. 



IV. AMPLIFICATION OF THE SUPERPOSED 
STATES VIA CONDITIONAL HOMODYNE 
DETECTION 

The fidelity between the state generated from our 
scheme and the ideal superposition of coherent states 
starts to decrease rapidly for a > 1. To realize surely 
macroscopic superposition or to apply these states to 
the coherent state superposition based quantum compu- 
tation scheme [17], the superposed states are required 
to have larger amplitudes. One approach for the pro- 
duction of such a state is to introduce PNRDs in our 
scheme (Fig. [3]) to generate a superposition of |^rn) 
and l^^+i) for large m. The other approach may be 
to apply the conditional amplification process proposed 
in [isl, which does not require PNRDs. It was shown 



\V2a) - I - V2a) 



that if one can prepare two inputs \a) + 



and 



they can be conditionally transformed to 
the state I7) + e^(^+^)| - 7), where 7 = ^/c^TW^ by 
using BSs, an auxiliary coherent state, and two on/off 
detectors. It could be used to amplify the initial state of 
5'(r)|l), which well approximates \C-{a)) for |ap < 1, 
and the scalability of the repetitive amplification process 
was discussed in detail in [20]. 

Application of our scenario to their scheme allows us 
to generate I7) + e*'^! — 7) with large 7 and arbitrary 
cp. Moreover, we will show in this section that the two 
on/off detectors used in the scheme of Ref. |0 can 
be simply replaced by a homodyne detector although the 
latter acts as a Gaussian operation. Note that it is not 
prohibited to transform non-Gaussian states to the other 
non- Gaussian states by only Gaussian operations. The 
conditional homodyne detection technique has been ap- 
plied to purify the coherent state superpositions ^ or 
the squeezed states suffering non-Gaussian noises (28|. 

The schematic of the conditional amplification via ho- 
modyne detection is shown in Fig.[5l[a). To explain how it 
works, let us see a simple example where we have 
and \C-{a)) as two input states. These states are com- 
bined by a balanced BS as 



{\a) + \-a)){\a)-\-a)) 

^ {\V2a) - I - V2a)m + \0){\V2a) 



(31) 
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FIG. 5: (a) A schematic of the conditional amplification of 
the coherent state superpositions with homodyne detection, 
(b) The three-step concatenation of the amplification process. 
The 8 initial states \a) + e^'^^ | — Q^) (i = I5 •••5 8) are prepared 
in certain phases Lp as illustrated in the left. 



where normalization factors are omitted for simplicity. 
Then one makes homodyne detection on one of the two 
modes. An ideal homodyne detection corresponds to a 
projection onto the quadrature eigenstate and for 
the measurement outcome x, one obtains the conditional 
output state 



(x|(|V2a) - I - 72a))|0) + (x|0)(|V2a) 
oc (e-(^-2")'/2 - e-(-+2-)'/2) |0) 

-"/^(\V2a)-\-V2a)). 



V2a)) 



(32) 



Here, conditioned on the outcome x = 0, the first term 
vanishes and one obtains the amplified state \^/2a) — \ — 
\/2a). More generally, the condition for the two inputs 
\a) + e^'^^'^l — a) to be amplified is (fi^(f>2 = tt. The 
amplified state has the phase ifi — if>2 which implies that 
one can choose arbitrary (p at the output. 
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For further amplification, the process should be con- 
catenated by carefully preparing the initial input states. 
Figure [Hb) is the schematic of the process to generate 
\2^/2a) — I — 2\/2a) by concatenating three amplification 
steps. The numbers represent the phase ^ of each state. 
The same rule can be applied in a straightforward way 
for the iterative generation of a superposition of large 
coherent states with arbitrary Lp. 

Homodyne detection is a well matured technique and 
very high quantum efficiency (j] > 0.99) has been 
achieved with current technology. It simplifies the exper- 
imental complexity, enhances the practical success prob- 
ability compared to a use of two imperfect on/off detec- 
tors, and thus will increase the total feasibility of the 
experimental demonstration. 



V. CONCLUSIONS 

In this paper, a novel scheme for the conditional gen- 
eration of a coherent state superposition with arbitrary 
ratio and phase has been proposed. The scheme uses 
a squeezed vacuum, beamsplitters, a coherent state an- 
cilla, and two on/off detectors. We first showed that 
c+|a) + c_| — a) for arbitrary {c+,c_} can be approx- 
imated by an appropriate superposition of single- and 



two-photon subtracted squeezed vacua with very high fi- 
delity {F > 0.99). 

Such a superposition of photon subtracted states is 
conditionally generated by using ideal photon number 
resolving detectors and a qubit ancilla bo\0) -\-bi\l). We 
have shown that this ideal scheme is also realized by 
the more practical scheme in which a displacement op- 
eration (coherent state ancilla) and on/off detectors are 
used. Even including realistic dark counts and low quan- 
tum efficiency (e.g. u = 10~^ and r] = 0.1), fidelities of 
more than 0.95 could be achieved with a highly trans- 
missive beamsplitter (T = 0.95) and the squeezing of 
r = 0.3 (~ 2.6 dB), which corresponds to the amplitude 
of a = 0.95 for the generated state. All these parameters 
are reasonably comparable with the recent experiments 
on single-photon subtraction 0, H, @ . 

We have also shown that our scheme is useful to 
simplify the setup of the conditional amplification of 
the coherent state superpositions originally proposed in 
[ii,[23. Our simplified version is quite feasible to demon- 
strate with current experimental techniques. Our scheme 
would also be useful to save the amount of required re- 
sources in other quantum information applications, in- 
cluding superposed coherent state based quantum com- 
puting pj,ll8|. 
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